METHOD FOR USING REMOTELY SENSED DATA 
TO PROVIDE AGRICULTURAL INFORMATION 

FIELD OF THE INVENTION 
5 The present invention is directed to a method for providing agriculture information on 

an area of interest, with the agricultural information being derived from the raw image data 
produced by a remote sensor associated with a remote sensing platform. 

BACKGROUND OF THE INVENTION 

10 The profit margins associated with the production and marketing of many crops have 

been generally decreasing for a number of years. This decrease in profit margins has caused 
many of the producers in the industry to turn to technology for information that facilitates the 
management of the crop producing land to increase the yield and/or reduce losses due to 
disease, drought and the like. 

15 Presently, one method for producing crop related information that allows the 

producers to manage crop producing land is to use remote sensing platforms (e.g., aircraft, 
satellites, and land vehicles) that are equipped with sensors that allow raw image data to be 
collected on the agricultural land of interest as the platform passes over the land. Typically, 
the raw image data is from one or more discrete bands in the electromagnetic spectrum that 

20 are capable of yielding agricultural information. For instance, an indication of the amount of 
chlorophyll in a plant may manifest itself in the "light" sensed from the land of interest in one 
or a combination of bands. This raw image data for the land of interest is subsequently 
processed to produce maps that provide the producers with information that can be utilized to 
manage crop production. For instance, a map that shows chlorophyll concentrations over the 

25 area of interest could be used by a producer to identify portions of the land that are under 
producing relative to other portions of the land and take appropriate action. 

SUMMARY OF THE INVENTION 
The present invention has recognized that much, if not all, of the agricultural 
30 information produced from the raw image data collected by remote sensing platforms is 
subject to sensor and/or atmospheric effects that lead to inconsistent data from one area of 
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interest to another area or interest, from one sensor to another sensor, and/or from one time to 
the next time. 

The present invention reduces the adverse effects of the sensor and/or atmosphere so 
that more consistent agricuhural information is produced thereby rendering the information 

5 easier for producers and others in the agricultural industry to utilize, to. other words, the 
invention produces more consistent agricultural information between different areas of 
interest. Consequently, if a producer has two different areas of land that are each planted 
with the same crop, a particular type of agricultural information produced by the invention on 
each of the two areas is substantially consistently represented, thereby facilitating the 

10 producers ability to compare the two crops to one another. 

In one embodiment, the present invention includes the steps of receiving information 
that defines the area of interest of the recipient of the agricultural information and receiving 
remote imaging data on the area of interest. The remote imaging data or raw data, which is 
typically produced by a sensor associated with a satellite or aircraft, is processed to produce 

15 reflectance factor data on the area of interest. To elaborate, the raw data is subject to effects 
that skew information on the area of interest. For example, the atmosphere between the area 
of interest and the sensor, whether satellite or aircraft based, if not taken into account, skews 
the agricultural information on the area of interest. By processing the raw data to produce 
reflectmce factor data, the effects of the atmosphere and/or other sources of skewing or error 

20 are significantly reduced. The method fiirther involves using the reflectance factor data to 
produce agricultural information on the area of interest. In one embodiment, the reflectance 
factor data is produced using agricultural scene objects. 

In one embodiment, the agricultural information is a green vegetation index (GVI) 
pixel map of the area of interest that provides an objective measure of the vegetation canopy. 

25 Each pixel on the map represents a particular portion of the area of interest. For instance, a 
pixel can represent a 5m by 5m area within the area of interest, to one embodiment, the pixel 
map utiUzes colors to indicate the GVI for the portion of the area of interest represented by 
each pixel. For example, brown can represent no vegetation, dark green can represent the 
densest vegetation canopy, and other colors can be used to represent canopies between these 

30 two extremes. 
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In another embodiment of the invention, the agricultural information is a pixel map 
that indicates changes in the GVI from date to date. Each pixel in the map is capable of 
representing a negative change, positive change or no change. The changes or lack of change 
associated with the land that each pixel represents are, in one embodiment, represented by 

5 colors. For example, shades of yellows and red can be used to represent negative changes in 
the GVI between two dates, shades of blue and green can be used to represent positive 
changes in GVI between two dates, and white can be used to represent no change in GVI 
between two dates. In one embodiment, the dates on which GVI information is produced 
that can be utilized to produce the map are within a week of each other. Consequently, the 

1 0 entity requesting the map is able to monitor changes between GVI maps on at least a weekly 
basis and potentially on a more frequent basis, such as day to day. 

In a further embodiment of the invention, the agricultural information is a vegetation 
index map other than a GVI map, such a map of the Difference Green Vegetation Index 
(DGVI), Perpendicular Green Vegetation Index (PGVI), Normalized Difference Vegetation 

15 Index (NDVI), Simple Ratio Vegetation Index (SRVI), Lifrared Percentage Vegetation Index 
(IRPVI), Soil Adjusted Vegetation Index (SAVI), Transformed Soil Adjusted Vegetation 
Index (TSAVI), Modified SAVI (MSAVI), or MSAVI2. Another embodiment apphes the 
reflectance factor data to produce two such maps that relate to different times and are used to 
generate a difference or change map. 

20 In yet a further embodiment of the invention, the agricultural information is a soil 

zone index (SZI) pixel map that characterizes an agricultural area of interest in terms of 
relative soil properties. In one embodiment, the relative soil properties that are characterized 
are a combination of surface soil wetness, soil texture, and organic matter content. The SZI 
map is considered most applicable to bare or mostly bare fields and is most likely to be 

25 utilized in making decisions regarding crop planning and management. 

In yet another embodiment of the invention, the agricultural information is a soil 
index map other than the SZI map, such as a map of the Soil Brightness Index. 

Another aspect of the invention involves the conveyance of the agricultural 
information produced from the reflectance factor data. To elaborate, the pixel map is 

30 typically produced and stored in a digital format that can be transmitted to the entity that 
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requested the agricultural information. In one embodiment, the map is transmitted to the 
requesting entity over the Internet. It is also feasible to record the map on a magnetic disk, 
CD, tape or other recording medium and then send the recorded medium to the requesting 
entity. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 illustrates a digital topographic map in which an area of interest (AOI) is 
located; 

Fig. 2 illustrate the digital topographic map of Fig. 1 with an overlying polygon that 
10 identifies the boundaries of the AOI within the area shown by the digital topographic map; 

Fig. 3 illustrates the digital topographic map of Fig. 1 with the overlying polygon that 
identifies the AOI, as displayed on a video monitor; 

Fig. 4 shows a panchromatic image of field an AOI that is a section of land with a 
number of fields; 

15 Fig. 5 illustrates the relationship of solar and atmospheric irradiances on the spectral 

radiance observed by a sensor directed at the earth's surface; and 

Fig. 6 illustrates a print out of a raster histogram tool that is capable of identifying a 
pixel in a radiance image that has a reflectance factor of effectively zero. 

20 DETAILED DESCRIPTION 

Generally, the present invention is directed to the production of agricultural 
information that is correlated to an area of interest (AOI), which is typically one or more 
fields. The process involved in producing the agricultural information includes receiving 
information that identifies the AOI, which is typically provided by the entity requesting the 

25 information. Also included in the process is receiving raw imaging data on the area of 
interest firom which agricultural information on the AOI can be derived. The raw data is 
processed to produce reflectance factor data, thereby substantially reducing the skewing 
effects attributable to the atmosphere and the sensor used to produce the data. The 
reflectance factor data is fiirther processed to produce agricultural information that is 

3 0 correlated to the AOI. 
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Having generally described the process for producing the agricultural information, an 
embodiment of the process is described in greater detail To provide agricultural infonnation 
on an AOI, the AOI is identified in a manner that allows the agricultural information that is 
to be subsequently produced to be correlated with the AOL One way to identify the AOI is 
5 to have the requester specify the boundaries of the AOI in terms of latitude and longitude. 
Specifying the AOI in this manner may be cumbersome for a requester. Consequently, to 
facilitate the identification of the AOI, the requester is provided with tools that facihtate the 
identification of the AOL Typically, the tools axe in the form of software that a requester can 
load onto a computer system. For AOIs in the United States, the tools include a digital 
10 topographic map product that is in the form of a digital raster of the USGS map(s) that 
M= encompass the AOL The USGS maps generally provide landmarks (roads, cities, waterways 

P- etc.) and township, range, and section information that is useful in identifying the AOI. An 

example of a topographic map product that encompasses an area of interest is shown in Fig. 
J: 1. It should be appreciated that the process is adaptable to maps other than USGS maps, 

15 provided the maps are capable of being correlated with the mapping information that is 
=^ associated with the raw data. 

The requester is also provided with an area polygon product for producing a vector 
Jl; data layer that allows the requester to delineate the AOI relative to the topographic map. 

D Stated differently, the area polygon product allows the requester to create a polygon that 

20 defmes the AOI relative to the topographic map by producing an overlay to the topographic 
map that delineates the AOL An example of a topographic map with overlying polygon that 
delineates the boundaries of an AOI is shown in Fig. 2. With respect to AOIs in the U.S., the 
requester is able to produce the polygon defining an AOI that is a section of land by 
specifying the township, range and section. The process is adaptable to other schemes for 
25 identifying parcels of land utilized in other countries. Figure 3 illustrates the digital 
topographic map of Fig. 1 with the overlying polygon that identifies the AOI, as displayed on 
a video monitor. 

Once an AOI is identified, the requester is further provided with a pan image product 
that is a panchromatic (black and white) image of the AOL This image is typically utiUzed 

30 by the requester to facilitate the identification of fields within an AOI in the maps that are 
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produced from the raw data and provide agricultural information on the AOL Figure 4 is an 
example of a panchromatic image that includes an AOI in the form of a section of land with a 
number of fields. 

A sensor, typically associated with a sateUite or aircraft, produces raw image data or 
5 raw image raster values (RVb) for each of the spectral bands that are needed to produce the 
agricultural information for the AOL More specifically, the sensor produces a raw image 
raster value for each of the electromagnetic bands of interest for each pixel in the scene (i.e., 
the area being viewed by the sensor) that includes the AOL In one embodiment, the sensor is 
capable of producing raw image raster values a blue hght band, green hght band, red Ught 
10 band and near infra-red band, which are also respectively referred to as band 1, band 2, band 
3 and band 4. It should also be appreciated that raw image raster values for a greater or 
lesser number of spectral bands may be required to produce a particular type of agricultural 
information. 

The raw image raster values for each band for each pixel that relates to the AOI is 
15 subjected to a sensor/radiance calibration step. For simphcity, the calibration is described 
with respect to the raw image raster value for a single pixel associated with a single band 
with the understanding that the production of certain agricultural information will require the 
processing of raw image raster values associates with multiple bands and multiple pixels. 
With reference to Fig. 5, the raster value associated with a pixel is represented by the 
20 following equation: 

RVb=Gb(x,y,t,Lb)Lb (1) 

where RVb is the raw image raster value for a band b for a pixel, Gb is the gain of the sensor 
25 for band b and varies over the images extent (x,y), over time (t), and perhaps with respect to 
the level of the spectral radiance being sensed (Lb). Gb also includes the effects of other 
sensor parameters, such as instantaneous field of view, and modulation transfer fimction. For 
simplicity, the ftinction Gb(x, y, t, Lb) is hereinafter referred to simply as Gb with the 
dependencies implied. Equation (1) can be rewritten as: 

30 
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Lb=RVb/Gb 



(2) 



Spectral radiance calibration of the RVb data is achieved by application of equation (2) to the 
RVb data, namely, dividing the RVb data by the value of Gb, which is known value for a 

5 given sensor. The appUcation of equation (2) to each pixel of an image for a particular 
spectral band produces a radiance image. 

The radiance image data is also commonly subject to positional or geometric skewing 
in terms of a given map coordinate system projection (mapping system and datum). Such 
skewing is or can be due to many factors (e.g., distortion in the camera system associated 

10 with the sensor, a tilt in the imaging system etc.). Correction involves the shifting of the 
pixels in the raster so that there is a substantially uniform distribution of pixels across the 
radiance image of the AOI in terms of the map coordinates of some selected mapping system 
and datum. If the data is subject to geometrical skewing, the data should be corrected so that 
the resulting map embodying the agricultural information conforms to the mapping system 

15 being utilized by the requester. Otherwise the resulting map may be very difficult for the 
requester to recognize as relating to the AOI. Such corrections to image data, such as 
radiance image data, are well-known in the art and, as a consequence, will not be elaborated 
upon further. 

After the raw image raster values (RVb) for the pixels in an image of the AOI and for 
20 each relevant band have been subjected to sensor/radiance calibration and/or correction for 
geometric skewing, atmospheric correction/reflectance calibration is performed on the 
radiance data Lb. As shown in Fig. 5, the total radiance encountered by the sensors is sum of 
two additive streams of radiant energy and can be represented as: 

25 Lb = Lpb + Lsb (3) 

where Lb is the total spectral radiance measured by the sensor, Lpb is the spectral radiance 
attributable to the scattering of electromagnetic radiation by gases and other materials in the 
atmosphere between the sensor and the surface of the earth; and Lsb is the spectral radiance 
30 attributable to diffuse reflectance by surface materials, such as soils and crops. Equation (3) 
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is rewritten to solve for the Lsb term, which relates to the soils and vegetation that are of 
interest, as follows: 

Lsb Lb ~ Lpb (^) 
Equation (3) can be represented by the following equation: 



where Eb is the spectral irradiance of the sun at the top of the atmosphere for band b; tsb is 
band transmittance (with the subscript s referring to the transmission loss between the sun 
and the surface of the earth); 6 is solar zenith angle; Eb,sky is the diffuse sky irradiance; Tob is 
the transmittance of the atmosphere between the surface of the earth and the sensor; and RFb 
is the reflectance factor for band b. Equation (5) applies to surfaces with Lambertian 
reflectance characteristics. Agricultural surfaces, however, have been found to exhibit non- 
Lambertian reflectance characteristics. This anisotropic behavior is addressed by replacing 
RFb with a bidirectional reflectance factor BRFb. Alternatively, RFb and n can be replaced 
with the bidirectional reflectance distribution factor BRDFb. The term RFb is retained 
hereinafter with the understanding that it refers to BRFb. To solve for RFb, equation (5) is 
rewritten as follows: 



Lb- 



(Eb tsb cos 6 + Eb,sky) tpb 



RFb + Lpb 



(5) 



7t 



RFb = 



7i(Lb-Lpb) 



(6) 



(Eb tsb cos 0 + Eb,sky) tob 



Equation (6) can be simplified to the following: 



RFb = (Lb - Lpb)mb 



(7A) 



where mb = 7r/(Eb tsb cos 0 + Eb,sky) tob 



(7B) 
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It is unlikely that all of the atmospheric parameters associated with mb and embodied in 
equation (7B) would be known as a sensor surveys a particular scene. However, the linear 
nature of equation (7 A) implies that only two independent items of information are needed 
about the surface reflectance within the scene to determine Lpb and mb. Specifically, if a dark 
5 object and a bright object can be found in the pixels of the radiance images with known 
reflectance factors, the two unknowns of equation (7 A), Lpb and mb, can be determined and 
the equation solved for RFb. To elaborate, if the known reflectance factor of the bright object 
is RFb,bnght and the known reflectance factor for the dark object is RFb,dark and the observed 
radiances for the bright object and dark object are respectively, Lb,bright and Lb,dark, then: 

10 

mb = (RFb,bright - RFb,dark) / (Lb,bnght - Lb,dark) (8) 
Lpb = Lb,dark - RFb,dark / mb (9) 

15 However, the estabUshment of objects within every pixel image that have reliable, known 
reflectance factors, i.e., the estabUshment of calibration targets in an AOI, is time consuming 
and expensive. 

The process utilizes the radiance image itself rather than any pre-established 
calibration targets in an AOI to determine mb and Lpb. To elaborate, with respect to 

20 determining Lpb, the radiance image produced according to equation (2) for an AOI almost 
always includes a pixel where the RFb is effectively zero or a feature with a known RFb at or 
about zero. Moreover, the smaller the resolution element (i.e., the smaller the area 
represented by a pixel), the greater the likelihood that there will be pixels with an RFb of 
effectively zero (shaded). The radiance image can be analyzed to identify a pixel or pixels 

25 that each have an RFb of effectively zero. Such analysis can be performed with a raster 
histogram tool in a GIS software package. Fig. 5 illustrates a radiance image with a pixel 
that has an Lb with a value of 296 and a histogram for the radiance image that indicates that 
the lowest Lb in the radiance image is 293. Consequently, the pixel that has an Lb with a 
value of 296 is a good candidate for representing Lpb within the image, i.e., the radiance 

30 attributable to the scattering of electromagnetic radiation by gases and other materials in the 
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atmosphere between the sensor and the surface of the earth. An example of a feature with a 
known RFb is non-turbid, deep open water, which has a known RFb of ne^ly zero in the near 
infra-red (NIR) band. 

For the pixel of the radiance image that is identified as having an RFb of 0 or 
5 approximately 0, equation (9) simphfies to: 

Lpb = Lb,dark (10) 

At this point, the value for Lpb has been determined. Further, the value of Lpb is applicable to 
the bands of interest. 
10 Rearranging equation (7 A) yields an expression for mb of: 

mb = RFb,bright /(Lb,bnght - Lpb) (1 1) 

Consequently, since Lpb is known from the radiance image (e.g. 296 in Fig. 5) and Lb,bnght 

15 will be known once a bright pixel is identified, mb can be determined if an object within the 
AOI can be identified with a known value for RFb,bnght. The process endeavors to identify for 
each band of interest an object that normally occurs in the agricultural AOIs that has a known 
or determinable RFb,bnght. Objects that meet this criteria may vary depending on the band of 
interest. For example, for the NIR band, dense herbaceous green vegetation (an agricultural 

20 scene-object) has a known reflectance factor of approximately 0.60. Consequently, for the 
NIR band, a pixel associated with dense herbaceous green vegetation is identified and the 
Lb,bright value for the pixel is used to determine m4 under equation (1 1). With respect to the 
blue, green and red Ught bands, as well as the NIR band, it has been found that dry, bare, 
agricultural soil, another agricultural scene-object, is a reasonably bright object. However, 

25 the RFb for dry, bare, agricultural soil is not fixed. In general, agricultural scene-objects are 
objects that relate either to the crop vegetation or the soil associated with crop field. 

To address the variable nature of the RFb for dry, bare, agricultural soil, a Line of 
Soils approach is employed. Specifically, there is a linear relationship between NIR 
reflectance factor of dry, bare, agricultural soils relative to red-light (RL) reflectance factor 

30 for the same soils that is expressed as: 
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RFnir,bas " Kir,rl + Snir,rl RFrl,bas 



(12) 



where the subscript has refers to bare agricultural soil, bn,r,ri is the intercept of the Line of 
5 Soils (i.e., the value of RFmr where RFh = 0 in the NIR versus RL feature space plot), and 
Snir,ri is the slope of the Line of Soils (in NIR versus Rl feature space). 

The Line of Soils in each of the three 2-space plots involving band 4 (i.e., the NIR 
band) are as follows: 



10 RF4,bas = b43 + S43 RFa.bas (13) 

RF4,bas = b42 + S42RF2.bas (1^) 

RF4.bas = b41 + S41 RFi,bas (15) 

15 



Once a pixel for bright, dry, bare, agricultural soil has been identified, the RF4,bas for the pixel 
is calculated using equation (7 A). The calculation of the RFs.bas utilizes the following known 
values of b43 and S43: 



20 b43 = 0.0; and (16) 

S43 = 1.1521 (17) 

These values may vary from one area to another. Further, these values vary based on the 
25 specific wavelengths and bandwidths of the NIR and RL bands employed. Consequently, 
these values are subject to modification. The values used for b42, S42, b4i and S41 are: 

b42 = 0.0 (18) 

30 S42 = 1.333 (19) 
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b4i = 0.0 



(20) 



S41 = 1.6667 (21) 

5 

These values are also subject to modification or refinement depending upon the situation. 
Based upon the substitution of the values for RF4,bas, b43 and S43 into equation (13), the value 
forRF3,basis: 

10 RF3,bas = RF4,bas/ 1.1521 (22) 



The value of m3 (i.e., the value of m for the red light band) is then calculated as: 

m3 = RF3,bas/(L3,bas - Lps) (2^) 

15 

Similarly, the values for ma and mi are calculated. At this point, the values for Lpb and mb in 
equation (7 A) have been determined for each band of interest. As a consequence, the values 
for RF for each pixel for each band in the radiance image of the AOI is determined to 
produce a reflectance factor image for each band. 

20 If the RF for dense herbaceous green vegetation in a particular scene is not 

approximately 0.60, errors in the estimates of m for the bands may result. In such a case, the 
calibration parameters Lpb and mb may be determined by other methods. One method, 
involves the use of non-agricultural objects in a scene m which the reflectance of the object is 
substantially constant, such as airport runways, a flat roof on a building and the like. If the 

25 scene has been previously rniaged, then a reflectance factor estimate exists that can be used 
in subsequent images to determine reflectance values if the RF for dense herbaceous green 
vegetation in the subsequent image is subject to question. Another method involves the use 
of calibration targets in the scene or other field measurements. Such an approach, in many 
cases, is not cost effective and care must be taken to use measurements at the same time that 
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the raw image data is produced and to account for solar illumination angles, observation 
angles etc. 

One or a combination of the reflectance factor images for the blue band (1), green 
band (2), red band (3) and NIR band (4) are used to produce an pixel map or image with 
5 agricultural information for the AOL One type of pixel image or map that is produced from 
the reflectance factor images is a calibrated green vegetation index (GVI) map that is 
indicative of the vegetation canopy in the AOL The equation for calculating the GVI based 
on RF values is: 

10 GVI = (RF4* - RF3*) (1 + L)F/(RF4* + RF3* + L) (24) 

where RF4* = RF4 - b43; (25) 
RF3* = S43 RF3; (26) 

15 

L = 0.5; and (27) 
F = 136. (28) 

20 The value of F is selected so that the value of GVI is about 100 when RF4 is 0.60 and RF3 is 
0.03. The parameter L is the soil adjustment coefficient. This results in GVI having a range 
of values from 0 to 100. A different range of values is also feasible by changing the value of 
F. The value of L, if the need arises, can also be changed. The values in the relevant GVI 
range are mapped to a color palette. In one embodiment, the lower values in the range are 

25 mapped to shades of browns, the upper values in the range are mapped to greens, and values 
in between the two extremes are mapped to other colors. It should be appreciated that the 
map can be a print out of the values for each pixel. However, such maps are typically more 
difficult to read than a map that maps the values to color. It should be appreciated that raw 
image data values or radiance values could be apphed to equation (24) if desired. 
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The GVI is believed to address certain deficiencies in the Normahzed Difference 
Vegetation Index (NDVI) and Soil Adjusted Vegetation Index (SAVI). Namely, under the 
NDVI and SAVI approaches the index value is only zero if RF4 is equal to RF3. This 
wrongly implies that the Line of Soils in RF4 versus RF3 space runs from (0,0) along a linear 
5 diagonal line of that space. In addition, NDVI is undefined when both RF4 and RF3 are equal 
to zero. NDVI is also highly sensitive to errors in RF4 or RF3 when the values are near zero. 
The NDVI also reaches a maximum value or saturates when the green vegetation reaches 
100% ground cover. As a consequence, the NDVI is largely incapable of giving any 
indication of the biomass that is produced after the vegetation reaches fiiU cover. The 
10 Transformed SAVI (TSAVI) approach accounts for the fact that the actual Line of Soils in 
RF4 versus RF3 space is given by a straight line that has a slope and intercept like that of 
P equation (13). The approach represented by equations (24), (25), (26), (27) and (28) is 

^; capable of accounting for the case in which the Line of Soils has a slope different than 1 and 

SI an intercept different than (0,0). 

Ill 15 Another type of pixel image or map that is or can be produced is a GVI change map 

:i - ii 

that reflects the changes in the GVI pixel images relating to two different times. Typically, 
O the requester desires the GVI change map to show the changes in GVI between two 

S5S23 

consecutive times when raw image data was obtained and from which a GVI map could be 
^ generated. However, if desired, a GVI change map showing the change in GVI between two, 

m 20 non-consecutive times can also be produced (i.e., when there are one or more times on which 
raw image data was obtained that are between the two times of interest). It should be 
appreciated that a GVI change map can be generated at any time provided that GVI data is 
available or can be produced for two different times. Typically, the GVI change map reflects 
the change in GVI between two times that are within one week of one another and more 
25 typically within three or four days of one another. In one embodiment, the raw image data 
for producing GVI images that are spaced four days from one another is accompUshed by 
the use of a constellation of satellites with the sateUites having orbit that result in one of the 
sateUite passing over a specific AOI every four days. If desired GVI change maps can be 
produced that use data for times that are more closely spaced in time. It should also be 
30 appreciated that images produced by other remote sensing platforms or combinations of such 
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platforms are also feasible. It has been found that producing GVI change maps using GVI 
data that is three to four days apart and providing the GVI change map to the requester within 
a few days of the end of the three or four day period provides the requester, especially when 
the requester is a farmer, with sufficient notice of changes in GVI to take appropriate action. 
5 In one embodiment, the GVI change data is mapped according to a scale that extends from 0 
to 200. A value of 100 indicates that there was no change in GVI between the two times that 
GVI pixel maps that are being used to generate the GVI change map were generated, A 
value of 200 reflects the maximum increase in GVI between the two times, and a value of 0 
reflects the maximum decrease in GVI between the two times. The scale is, in tum mapped 

10 into a color palette. In one embodiment, GVI change values between 0 and 94 are mapped to 
shades of red and yellow; a GVI change value of 95-105 is mapped to white; and GVI 
change values between 106 and 200 are mapped to shades of blue and green. 

Another type of pixel image or map that is or can be produced from the reflectance 
factor images is a soil zone index (SZI) map that is indicative of differences in soil 

15 properties. If the GVI pixel map is generated, the equation for determining the SZI is: 

SZI = 250 [ sqrt( (RF3* + L/2f + (RF4* + L/2f ) - 0.35 ] (29) 

In practice, the range of SZI extends from 0 to 100. Low values of SZI (0-16) are associated 
20 with dark areas that are in shadows; values of SZI near 17 are associated with wet soils, silty 
loam soils, and/or soils having high organic matters content; and values of SZI near 100 are 
indicative of dry soils with sandy textures and low in organic content. The SZI values are 
mapped to a color palette to produce a colorized SZI pixel map. A different scale is also 
feasible and can be reahzed by changing the multiplication constant in equation (29). It 
25 should be appreciated that the map can be a print out of the values for each pixel. However, 
such maps are typically more difficult to read than a map that maps the values to color. It 
should also be appreciated that raw image data values or radiance values can be applied to 
equation (29) if desired. 

The reflectance factor data, while discussed hereinabove with respect to the 
30 production of a GVI map, is also capable of being apphed to known vegetation indexes, such 
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as the Difference Green Vegetation Index (DGVI), Perpendicular Green Vegetation Index 
(PGVI), Normalized Difference Vegetation Index (NDVI), Simple Ratio Vegetation Index 
(SRVI), Infrared Percentage Vegetation Index (IRPVI), Soil Adjusted Vegetation Index 
(SAVI), Transformed Soil Adjusted Vegetation Index (TSAVI), Modified SAVI (MASVI), 
5 and MSAVI2. Moreover, the reflectance factor data is capable of being used in the 
production of change maps for each of the indexes. 

The reflectance factor date, while discussed hereinabove with respect to the 
production of an SZI map, is also capable of being appHed to known soil indexes, such as the 
Soil Brightness Index (SBI). 
1 0 Whatever agricultviral information map is produced using the reflectance factor values 

is conveyed to the requester. In one embodiment, the map is transmitted to the requester by 
conveying the map over the Intemet. Typically, the map is conveyed in a compressed 
format. Once received, the requester is able to produce an layered image of the AOI with the 
bottom most layer being the digital topographic map in which the AOI is located, the next 
15 layer bemg the vector polygon that defines the boundaries of the AOI, and the final layer 
being the agricultural information map that overiays the polygon defming the AOI. It is also 
possible to convey the agricultural information map to the requester in other ways. For 
instance, the agricultural information map can be recorded on a magnetic disk, CD, tape or 
other recording medium and mailed to the requester. If needed the recording medium can 
20 also include the digital topographic map and the vector polygon layers also. It is also 
possible to simply produce a hard copy of a composite map in which tiie agricultural 
information map overlays the relevant topographic map and vector polygon that defines the 
AOI and then mail the hardcopy the requester. The hard copy can also be faxed or otherwise 
electronically sent to the requester. 
25 The embodiment of the invention described hereinabove is intended to explain the 

best mode known of practicing the invention and to enable others skilled in the art to utiUze 
the invention. 
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